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Background 

When a single alloy is present in a 
flame, the identity of its elemental con- 
stituents, and therefore the alloy 
itself, can readily be determined using 
spectroscopic techniques. If multiple 
alloys are present in the flame, the 
identification process is again rela- 
tively simple if the alloys are composed 
of different elements. However, many 
alloys contain some of the same elemental 
species, their differences being various 
weight fractions of the elements. If 
numerous alloys are present in a flame, 
and these alloys contain some of the same 
elemental constituents, the spectra of 
the individual elements are summed, and a 
direct determination of identity and 
amount of the alloy cannot be made solely 
with spectroscopic techniques. 


Abstract 

The mathematical derivation for a 
method to determine the identities and 
amounts of alloys present in a flame 
where numerous alloys may be present is 
described. This method is applicable if 
the total number of elemental species 
from all alloys that may be in the flame 
is greater than or equal to the total 
number of alloys. Arranging the atomic 
spectral line emission equations for the 
elemental species as a series of simulta- 
neous equations enables solution for 
identity and amount of the alloy present 
in a flame. This technique is intended 
for identification and quantification of 
alloy content in the plume of a rocket 
engine. Spectroscopic measurements 
reveal the atomic species entrained in 
the plume. Identification of eroding 
alloys may lead to the identification of 
the eroding component. 

Nomenclature 

A transition probability 

Av Avogadro's number 

E energy 

g statistical weight 

h Planck's constant 

k Boltzmann's constant 

mw molecular weight 

$ radiant power 

T temperature 

v frequency 


The emission spectrum of an alloy is 
the sum of the line emissions from its 
elemental constituents. Identifying key 
spectral lines from the various atomic 
species enables identification of the 
elements. Since the radiance of a spec- 
tral line is related to the number of 
atoms in the emitting system, atomic con- 
centration can also be determined. If 
numerous alloys contain some of the same 
elements, the atom population for the 
various elements may be due to contribu- 
tions of the element from numerous 
alloys. The resulting spectral radiance 
that is measured from an atomic species 


may be due to contributions of the ele- 
ment from numerous alloys. 

Emission spectroscopy of a rocket 
engine exhaust plume is being investigat- 
ed to determine if engine degradation can 
be monitored from measurements of mate- 
rial erosion. 1 ' 4 Knowledge of material 
erosion from internal components may pro- 
vide insight into the engine's current 
health and remaining life. In the space 
shuttle main engines (SSME) , it has been 
determined that material does in fact 
erode from internal engine components. 
Eroded particles are carried along with 
the hydrogen and oxygen propellants into 
the combustion chamber, where they are 
burned (dissociated to atomic species), 
thermally excited, and emit their charac- 
teristic spectrum. Figure 1 displays the 
propellant flow path within the SSME. 

Note that the propellants travel through- 
out the entire engine, with the ultimate 
destination being the combustion chamber. 
Figure 2 displays a typical SSME emission 
spectrum. Atomic lines can readily be 
identified in this figure. Alloys from 
which the SSME is fabricated are listed 
in Table 1. This table lists the alloys, 
their elemental constituents, and the 
weight fraction of these elemental con- 
stituents. 5 As can be seen from this 
table, there are numerous alloys in the 
engine, many of which contain some of the 
same elements. Table 2 is a list of the 
SSME components in the propellant flow 
path and the alloy from which they are 
fabricated. Note that the propellants 
and impurities in the propellants are not 
included in the above tables. 

This paper presents a technique that 
enables determination of the identities 
and amounts of alloys that are present in 
a flame. The mathematical derivation, 
and errors associated with this technique 
are discussed. Note that the primary 
purpose of this paper is to document the 
mathematics describing this technique, 
with minimal consideration given to 
implementation, calibration, and opera- 
tion of this capability. 


constituents, weight fraction of their 
elemental constituents, spectral charac- 
teristics of these elements, and the tem- 
perature at the location of the emitting 
system, a series of simultaneous equa- 
tions may be developed and solved. The 
solution to these equations yields the 
amount and identity of the alloy in the 
flame. Current state of the art in 
rocket engine diagnostics using plume 
spectral measurements enables identifica- 
tion and quantification of the atomic 
species. Application of the technique 
presented in this paper to a rocket 
engine enables identification and quant- 
ification of eroding alloys, which pro- 
vides information that can lead to the 
identification of the eroding engine com- 
ponent. Note that this technique is not 
specific to rocket engine diagnostics; 
it may be used in other applications 
requiring alloy identification and 
quantification. 


Theory 

Radiance from atomic line emissions 
is a function of several variables. In 
this paper, the general expression for 
line radiance as a function of atomic 
population densities shall be considered. 
It is assumed that the collection of 
atoms in the emitting system are in ther- 
mal equilibrium, and that there is a low 
population density of the emitting atoms. 
It is also assumed that the radiance of 
the line is sufficient to enable spectro- 
scopic measurements. The radiant power 
of an atomic line is given by the equa- 
tion: 6 


$ = 


A j>ji9j ne ~ Ej/kT 

zvn 


( 1 ) 


where ji represents one of the allowed 
transitions for an atomic species from an 
upper energy state j to a lower energy 
state i, and n is the total number of 
free atoms in the emitting system. 
Selection of a particular transition 
defines the frequency i/ and therefore 
the wavelength of the emitted photons. 


From knowledge of which alloys may 
be present in the flame, their elemental 
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Numerical values for the transition prob- 
ability, statistical weight and upper 
energy state for the selected transition 
can be found in reference tables. The 
partition function can either be found in 
reference tables 8 or calculated from the 
following equation: 


molecular weight and multiplying by 
Avogadro's number. Therefore, Eq. (5) 
can be rewritten as: 


x a Av 





mw,. 


( 6 ) 


Z(t) ■ Jj, gj.-W" (2) 

The solution to Eq. (1) yields the 
radiation at a given wavelength from a 
collection of atoms at a particular tem- 
perature. Utilizing emission spectro- 
scopy techniques, this radiance can be 
measured. Therefore, at a known tempera- 
ture for an allowed transition, Eq. (1) 
can be solved for the total number of 
free atoms in the emitting system. The 
number of atoms for the various elemental 
species can be used to determine the 
amount of the alloy that is present in 
the flame. To simplify the writing of 
Eq. (1), the following definition will be 
made: 

n A ji h2/ ji9j e " Ej/kT (3) 

= — zrn — 


where N a is the total number of atoms 
in the amount x of alloy a. Writing 
similar equations for m different 
al loys yields: 


Nj = XjAv 


N 2 = x 2 Av 


N m = x m Av 


w 12 

+ • • • 

+ V 

mw 2 

mw p 

w 22 
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+ w 2p 

mw 2 

mw p 
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(7-1) 


(7-2) 


(7-m) 


where R represents the radiant power 
from a single atom. Substituting Eq. (3) 
into Eq. (1) yields: 

$ = Rn ( 4 ) 

Consider an unknown amount x, of 
alloy a. The alloy consists of a ele- 
mental constituents, with the weight 
fractions of these elements equal to w p 
w 2 , . . ., w . It is assumed that the 
exact composition and weight fractions 
are known. The total weight w of alloy 
a can be expressed as: 

w a = x a (w 1 + w 2 + . . . + w,) (5) 

This equation can be expanded to 
calculate the total number of atoms in 
the alloy based on the total number of 
atoms of each elemental species by divid- 
ing the weight of an element by its 


where N : through N are the number of 

atoms in alloys 1 through m; x, 
through x m are the amounts of alloys 1 
through m; mw, through mw are the 
molecular weights of elemental species 1 


through p, and w n through vil are 
the weight percentages of the elemental 


species in the alloys. Note that the 
number of elemental species in Eqs. (7-1) 
to (7-m) has changed from a to p. 

While c represents the total number of 
elemental species found in a single 
alloy, p represents the total number of 
different elemental species that can be 
found in all alloys 1 through m. With 
this type of notation, the weight per- 
centage of an element in the alloy is set 
equal to zero if the particular element 
is not present in an alloy. 


By summing vertically down the col- 
umns of the above equations, the total 


3 


t 


number of atoms of each elemental species 
can be calculated: 

n l = + X 2 W 21 + ••• + x m w ml) (8 - l) 

UIW| ' ' 


n 2 = ( x l w 12 + X 2 W 22 + • • • + x m w m2) (8"2) 

IIIW2 ' 


n P = ^ r ( X 1 W 1 p + X 2 w 2p + • • • + x m w m p) (8-p) 

lilWp 


where n, through n are the total 
number of atoms of the each elemental 
species contributed from alloys 1 
through m. Substituting these values 
into Eq. (4) and calculating a value for 
R from Eq. (3) enables radiance calcula- 
tions of the p elemental species as a 
function of the amount x of the alloy. 
With this substitution, the following 
series of equations can be written: 


RjAv 


*1 = “iiSr( XlWl1 + X 2 W 2! + •” + x m w ml) ( 9_1 ) 


mwj 


$0 = 


R 2 Av 

mw 2 


( X 1 W 12 + X 2 W 22 + . . . + x m w m2 ) (9-2) 


R Av 

$p = li^( XlWl P + X2W2p + ••• + XmWmp ) (9 " p) 

At this point, consideration must be 
given to the known and unknown variables 
in Eqs. (9-1) to (9-p). The radiant 
power $ is measured spectroscopically 
from the flame, and therefore is known. 

It is assumed that the temperature of the 
emitting system is known, resulting in 
known values for Rj through R as 
defined by Eq. (3). The weight percent- 


ages of the elemental species in the var- 
ious alloys are obtained from metallurgi- 
cal handbooks, while the elemental molec- 
ular weights and Avogadro's number are 
known. The unknowns are the amounts of 
the alloys x, through x . Therefore, 
examination of Eqs. (9-1) to (9-p) indi- 
cates that there are p equations with 
m unknowns. If p is greater than or 
equal to m, a sufficient number of equa- 
tions are available to solve for all 
unknowns. If m is greater than p, the 
equations cannot be solved. 


Mathematical Implementation 


Assuming that there are a sufficient 
number of equations to solve for the 
unknowns, (the total number of different 
elemental species is equal to or greater 
than the number of alloys) the following 
matrix expression of Eqs. (9-1) to (9-p) 
can be written: 


R l w ll 

R 1 w 21 

mwj 

m 1 

R 2 w 12 

R 2 W 22 

mw 2 

mw 2 

R l w lp 

R l w 2p 

mwp 

mw p 



The above matrices can be expressed 
as [R] [X] = [$] , where the R matrix 
contains the R^ w l4 Av/mw 1+ terms, the 
X matrix contains tne x 1 through x m 
terms, and the $ matrix contains the 
through $ terms. Values for the 
R matrix are derived by simple calcula- 
tions, values for the $ matrix are 
measured spectroscopically from the 
flame, and the X matrix contains the 
unknown mass and by definition the iden- 
tity of a particular alloy. If the 
total number of elemental species is 
greater than the number of alloys, (p-m) 
equations can be eliminated. After 
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elimination of the excess equations, the 
R matrix is square and can be inverted. 
The solution to the X matrix is 
[X] = [R]' 1 [$]. The matrix inversion and 
subsequent matrix multiplication can be 
performed on a computer. Given that the 
[R] matrix can be inverted, a solution 
will be produced. The accuracy of the 
solution is limited by the accuracy of 
the values in the [R] and [$] 
matrices . 

Examination of Table 1 reveals that 
there are 17 alloys consisting of 20 ele- 
mental species in the SSME. Therefore, 
three elemental species do not need to be 
spectroscopically monitored for alloy 
identification and quantification. The 
selection of which elements to eliminate 
is based on the relative abundance and 
emission intensity of the species. Since 
zirconium, lanthanum, and nitrogen appear 
in trace quantities in a very limited 
number of alloys and are relatively weak 
emitters, 9 they should be the species 
that are eliminated. 


Discussion 

Although the mathematics described 
above are relatively straightforward, 
implementation of this technique requires 
additional considerations . These consid- 
erations include the errors that are 
inherently associated with the technique, 
and the sensitivity, dynamic range, and 
calibration of the instrumentation used 
to measure the emission spectrum. Fol- 
lowing is a discussion of the errors and 
instrumentation issues. While there may 
be numerous sources of error, the errors 
do not appear to be insurmountable indi- 
vidually or collectively. 

The derivation of this technique is 
based on spectroscopic, metallurgical, 
and mathematical considerations. Errors 
associated with the implementation of 
this technique to the analysis of a rock- 
et engine exhaust plume are the result of 
the limitations and assumptions in spec- 
troscopic theory and spectroscopic meas- 
urement techniques, uncertainties in 


alloy content, and the assumption that 
all constituents of the alloy wear at a 
constant rate. Spectroscopically, the 
assumptions that may lead to error 
include: (1) the temperature is known 

and uniform throughout the emitting sys- 
tem, (2) the radiance can be accurately 
measured regardless of the number of 
atoms and radiance from the emitting sys- 
tem, and (3) that exact values for the 
transition probabilities are known. Met- 
al lurgical ly, the assumption that the 
exact composition and weight percentages 
of the elements in the alloys are known 
may lead to error. The assumptions that 
the various elemental species in the 
alloy erode evenly and at the same rate 
is also a potential source of error. 
Following is a discussion of the various 
errors that may be caused by the 
assumptions. 

Temperature dependence of the radi- 
ance from atomic emissions can be seen in 
Eq. (1). For the SSME, temperatures in 
the combustion chamber and Mach disk 
regions are well established from analyt- 
ical calculations. 10,11 The assumption 
that the emitting atoms are in thermal 
equilibrium is valid since the tempera- 
tures in the chamber or the Mach disk are 
reasonably uniform. Temperature varia- 
tions occur near the walls of the cham- 
ber, and near the edges of the Mach 
structure. If spectral emissions from 
the center of the combustion chamber or 
center of the Mach structure are measured 
and used in the Eq. (9), thermal equilib- 
rium and uniform temperature is assured. 
Note that using a limited field of view 
introduces the assumption that metallic 
species entrained in the plume are homo- 
geneously mixed within it; and that 
observation of a known percentage of the 
plume accurately reflects the overall 
composition of the entire plume. Meas- 
urements from a limited field of view 
requires that a scaling factor be incor- 
porated into Eq. (9) to account for the 
hot gases and therefore spectral emis- 
sions that are not being observed. 

To gain insight into the temperature 
effects with respect to the SSME, the 
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Boltzmann distribution (temperature 
dependent terms in Eq. (1)) were evalu- 
ated at 3000 and 3149 K for chromium and 
titanium. Numerical values for the 
energy levels, statistical weights, and 
partition functions were obtained from 
Refs. 7 and 8 respectively. It was cal- 
culated that for the 5 percent tempera- 
ture change, the radiant power increased 
46.0 percent for chromium, and 
83.6 percent for titanium. As seen from 
Eq. (1), the temperature term appears in 
an exponential, and the significance of 
this term depends on the slope of the 
exponential at the temperature regime of 
the emitting system. For the SSME tem- 
perature regime, temperature is critical 
for accurate measurements, and a 
5 percent temperature deviation would 
yield unusable results for this 
technique. 

It is assumed that an adequate num- 
ber of free atoms of the species being 
spectrally monitored are available within 
the observed region to enable measure- 
ment. A single or very few emitting 
atoms of a particular species will not be 
seen in an environment consisting of many 
free atoms from numerous species due to 
either the low radiance from the few 
atoms, poor sensitivity of the detector, 
or low transition probability of the 
transition that is being monitored. If 
an extremely large number of atoms of a 
particular species are present in the 
flame, the flame is optically thick 
rather than optically thin, requiring 
that the self absorption terms be 
included in Eq. (1). Therefore, the 
population of the metallic species must 
be sufficient to provide radiance in 
excess of that of the background chemi- 
luminescence from the combustion process, 
yet not great enough to cause self 
absorption. In the SSME, the population 
of the eroded metallic species from the 
alloys is much less than the population 
of the oxygen and hydrogen propellants 
(unless. a major failure is occurring). 

Error may result from the transition 
probabilities and their associated uncer- 
tainties. Radiance emitted from an 
excited atom is directly proportional to 


the transition probability of the speci- 
fic transition that is being monitored. 
Low transition probability could lead to 
error due to the low overall radiance 
from the atoms in the emitting system. 

The uncertainty associated with the vari- 
ous transition probabilities may also 
cause errors. Typically, values for the 
transition probabilities are listed with 
an associated uncertainty band. These 
uncertainties range from 1 percent up to 
50 percent for the various transitions as 
seen in Ref. 7. Selection of a transi- 
tion with a low uncertainty will minimize 
error; however this may be difficult, 
since not all species have been ade- 
quately studied to generate high accuracy 
transition probabilities. 

Sensitivity, signal to noise ratio, 
dynamic range, and resolution of the 
spectrometer may limit the accuracy of 
this technique. The sensitivity and sig- 
nal to noise ratio dictate the minimum 
radiance level that can be measured. The 
dynamic range limits the maximum radiance 
level that can be measured, assuming the 
spectrometer is optimized for low radi- 
ance levels. Spectrometer resolution 
must be adequate to separate closely 
spaced spectral lines. It was assumed 
that spectral emissions from all elemen- 
tal species that may be in the plume can 
be measured. This assumption is correct 
if the detector is sensitive enough to 
measure potentially low radiance levels. 
From Table 1 it can be seen that some 
alloys contain trace amounts of certain 
elements. It may be extremely difficult 
to measure spectral emissions from these 
trace amounts unless a large quantity of 
the alloy is eroding, or the transition 
that is being monitored has a large tran- 
sition probability. Past work 4 has indi- 
cated that large amounts of some elements 
must be present in the plume for observa- 
tion due to low transition probabilities. 
The worst case scenario is the low weight 
percentage species that has low transi- 
tion probabilities. The lack of sensi- 
tivity may be overcome with higher 
sensitivity detectors, or by using atomic 
absorption rather than atomic emission 
spectroscopy. 
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Spectrometer calibration must also 
be considered since absolute, rather than 
relative spectral emissions are being 
measured. Calibration must take into 
account spectrometer throughput, effi- 
ciency of the detector at all observed 
wavelengths, and view angle of the plume. 
It may be possible to calibrate the spec- 
trometer at a few spectral lines using 
calibrated lamps and use these measure- 
ments as the basis for calibration of the 
spectrometer. 

The assumptions that the composition 
of the alloy is known exactly and that 
the elemental constituents of the alloy 
are homogeneously mixed may lead to 
error. Typically, the composition of an 
alloy is expressed as the nominal, rather 
than exact, weight percentages of its 
elemental constituents. The exact per- 
centages of the elements are not known 
unless a thorough analysis of the alloy 
is made. Also, the alloys may contain 
trace impurities. If an element is a 
trace impurity in one alloy but a consti- 
tuent in another alloy, errors may 
result. For critical applications such 
as the SSME, it is presumed that the 
exact composition of the alloys are 
known. 

To enable a mathematical solution, 
the number of elemental species must be 
greater than or equal to the number of 
alloys that are found in the engine. 
Examination of Table 1 will verify that 
this is the case for the SSME. Note that 
an appropriate set of equations must be 
generated for the temperature of the 
emitting system. 


Concluding Remarks 

The mathematical derivation of a 
technique that enables determination of 
the identities and amounts of alloys 
present in a flame where numerous alloys 
may be present has been described. This 
technique is particularly useful when the 
same elemental species are found in many 
of the different alloys. The intended 
application of this technique is to 
enable identification and quantification 


of eroding material from internal compo- 
nents of a rocket engine based on its 
plume spectral emissions. Identification 
of the eroding alloy may lead to identi- 
fication of the eroding component, and 
quantification of erosion may lead to 
determination of the remaining life of a 
component. 

Requirements for implementation of 
this technique include knowledge of which 
alloys may be present in the flame, their 
elemental constituents, weight fraction 
of the elemental constituents, spectral 
characteristics of these elements, tem- 
perature at the location of the emitting 
system, and measured radiance values from 
the atomic emissions in the emitting sys- 
tem. Arranging the atomic emission equa- 
tions for all the elemental species in 
the alloys as a series of simultaneous 
equations enables solution for the iden- 
tity and amount of the alloy in the 
flame. To enable a mathematical solu- 
tion, the total number of elemental spe- 
cies from all alloys must be greater than 
or equal to the number of alloys that are 
found in the engine. An appropriate set 
of equations must be generated for the 
temperature of the emitting system, and 
that if the temperature changes, a new 
set of equations must be generated. Note 
that an "alloy" may consist of a single 
element . 

The viability of this technique is 
limited by the inherent errors and 
assumptions that went into its deriva- 
tion. Although there are numerous 
sources of errors, they do not appear to 
be insurmountable. Critical sources that 
may lead to the most significant error 
are the uncertainty and likely 
fluctuation of the temperature in the 
emitting region, and the uncertainty 
associated with the transition 
probabi 1 ities. 

It is assumed that an adequate num- 
ber of free atoms of the species being 
spectrally monitored are available within 
the observed region to enable measure- 
ment, and that spectral emissions from 
all elemental species that may be in the 
plume can be measured. The population of 
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the metallic species must be sufficient 
to provide radiance in excess of that of 
the background chemiluminescence from the 

combustion process, yet not great enough 
to cause self absorption. The sensitiv- 
ity, signal to noise ratio, dynamic 
range, and resolution of the spectrometer 
that is measuring the spectral emissions 
must be adequate for the widths and peak 
intensities of the spectral lines being 
measured. 

The composition of the alloys must 
be known, and the elemental constituents 
of the alloy must be homogeneously mixed. 
The assumption that the exact composition 
of the alloy is known may lead to error, 
since alloy composition is typically 
expressed as the nominal weight percent- 
ages of its elemental constituents. 

Also, the alloys may contain trace impu- 
rities. If an element is a trace impuri- 
ty in one alloy but a constituent in 
another alloy, errors may result. The 
elemental constituents of the alloy must 
be homogeneously mixed, or the erosion 
rate of the elementals will not be 
constant. 
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TABLE 1. - NOMINAL ELEMENTAL COMPOSITION OF SSME ALLOYS IN WEIGHT PERCENTAGES 
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TABLE 2. - SSME COMPONENTS AND MATERIALS (PARTIAL LIST) 


Materi al 

Found in SSME part 

A-286 

Fuel and oxidizer preburner 

Main fuel, oxidizer and coolant valve 

AISI 440A 

Fuel, oxidizer and coolant valve body 

AISI 440C 

HPOTP and HPFTP bearings 

Armco 21-6-9 

Oxidizer pump discharge flex joint 
LPFP discharge flex joint 
Fuel bleed line flex joint 

304 CRES 

Main injector heat shield 

316L CRES 

Heat exchanger 

347 CRES 

Main injector faceplate 

Haynes 188 

LOX posts 

Heat exchanger 

Fuel bleed line flex joint 

Incoloy 903 

External gimbal ring flex joint 
Hot gas manifold liner 
Heat exchanger 
HPOTP primary turbine seal 
HPOTP secondary turbine seal 
HPOTP turbine interstage seal 

Inconel 625 

Temperature sensors 
Fuel preburner faceplate 
Oxidizer preburner faceplate 
Heat exchanger 
Hot gas manifold 
Heat exchanger 
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TABLE 2. - CONCLUDED. 


Material 

Found in SSME part 

Inconel 718 

Oxidizer pump discharge flex joint 

LPFP discharge flex joint 

Fuel bleed line flex joint 

HPOTP intermediate seal 

HPOTP boost pump labyrinths 

HPOTP pump impeller 

HPOTP pump stator 

Main combustion chamber shell 

External gimbal ring flex joint 

Oxidizer tank pressurization duct flex joint 

Hot gas manifold 

Heat exchanger 

K-monel 

HPOTP primary oxidizer seal 

MAR-M246+Hf 

HPOTP turbine blades 
HPFTP turbine blades 
HPOTP turbine tip seals 

NARloy A 

Preburner baffles 

NARloy Z 

Main combustion chamber 

Rene 41 

HPOTP turbine tip seal 

Stel 1 ite 
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Figure 1 . — SSME propellant flow paths. 
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Figure 2 — Typical SSME emission spectrum. 
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